Class of consistent fundamental-measure free energies for hard-sphere mixtures 
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In fundamental-measure theories the bulk excess free-energy density of a hard-sphere fluid mixture 
is assumed to depend on the partial number densities {pi\ only through the four scaled-particle- 
theory variables {£ a }, i.e., &({pi}) — > ^({Cq})- By imposing consistency conditions, it is proven here 
that such a dependence must necessarily have the form = — fo fn(l — &) + i &(y)tii&/{l — £3), 

where y = £|/127r£i(l — £3) is a scaled variable and ^(y) is an arbitrary dimensionless scaling 
function which can be determined from the free-energy density of the one-component system. Ex- 
tension to the inhomogeneous case is achieved by standard replacements of the variables {£, a } by 
the fundamental- measure (scalar, vector, and tensor) weighted densities {n a (r)}. Comparison with 
computer simulations shows the superiority of this bulk free energy over the White Bear one. 

PACS numbers: 05.70. Ce, 61.20.Gy, 65.20. Jk 



Introduction. Rosenfeld's fundamental-measure the- 
ory (FMT) 0, 0, together with its variants is 
perhaps the most successful density-functional approach 
for the description of inhomogeneous hard-sphere (HS) 
systems. Recent reviews can be found in Refs. (Tll - 
Il3j | - FMT is based on the ansatz that the excess free- 
energy density ^(r), which in general is a functional of 
the partial number densities Pi(r'), can be approximated 
by a function of a reduced number of weighted densities 
G3 K(r)}, i.e., $[r,{ Pi (r')}] -> $({n a (r)}). In order 
to choose a specific form for $({n a (r)}), the bulk free- 
energy density $({pi}) —> &({£,a}), where {£ a } are the 
so-called scaled-particle theory (SPT) variables, plays a 
crucial role. In Rosenfeld's original formulation [Tl|, 
application of two thermodynamic consistency conditions 
[cf. Eqs. © and © below] yields the SPT free energy. 
Later on, Roth et al. [|[ sacrificed consistency in favor of 
accuracy and proposed the so-called White Bear (WB) 
version, which is consistent with the successful Boublik- 
Mansoori-Carnahan-Starling-Leland (BMCSL) equation 
of state (EOS) pHH via Eq. © but not via Eq. ©. 
More recently, Hansen-Goos and Roth [h], O, E3 have 
proposed a modification of the WB functional (the so- 
called WB mark II) such that the inconsistency with Eq. 
([5]) disappears in the one-component limit. 

The main aim of this Rapid Communication is to show 
that, if the consistency conditions © and are com- 
plemented with a recently derived scaling relation fl8| . 
the bulk free-energy density of the mixture is unambigu- 
ously expressed in terms of the bulk free energy of the 
one- component system. If, in particular, the Carnahan- 
Starling (CS) expression [19( is chosen as input, the out- 
put results for the bulk mixture are more accurate than 
those resulting from the WB free energy. 

Bulk free-energy density. Let us consider a three- 
dimensional (additive) HS fluid mixture with (bulk) par- 
tial number densities pi and (bulk) total number density 
p = J2iPi- The moments of the size distribution arc 



M a = ^2iXiaf, where Xi = pt/p and Oi are the mole 
fraction and diameter, respectively, of spheres of compo- 
nent i. 

Let ( { Pi } ) be the bulk excess free-energy density, in 
units of fc^T, where ks is the Boltzmann constant and 
T is the absolute temperature. The free energy is said 
to be "truncatable" [20l - [22| if $({pi}) depends on {p t } 
only through the total packing fraction 77 = (7r/6)pM 3 
and a finite number K of moments {Mi, M2, . . . , Mk}- 
According to FMT arguments, the relevant number of 
moments for three-dimensional systems is K = 3. In 
that case, instead of the set of variables {77; Mi, M2, M3}, 
one can alternatively take the four SPT variables £0 = P, 
fi = \pM u ^ = 7rpM 2 , and £ 3 = f pM 3 . 

By requiring consistency with the limit where one of 
the components is made of point particles, it has been 
proven in Ref. [IH that a truncatable free energy (with 
K = 3) must necessarily have the scaling property 



$({&»}) = -£oln(l - 6) + Zi&H(ta,z), 



(1) 



where we have called z = £f /127r£i and the dimensionless 
scaling function H(^,z) remains undetermined. Equa- 
tion {l} imposes an important constraint on the func- 
tional dependence of $ on the four SPT variables 
since the unknown function H depends on two variables 
only [l8j . The main goal of this Rapid Communication is 
to prove that an extra consistency condition further con- 
strains $ to an unknown function of one variable only. 

First, we consider the standard thermodynamic rela- 
tion 



/9 + /T 



d $ 
dp p 



£0 - * + J2 & 



(2) 



where p is the pressure and /3 
SPT condition is d H 



1/ksT. An independent 



= lim 



Pp7 

™f/6 



(3) 



' andrcs@uncx.es http://www.uncx.es/cwcb/fistcor/andrcs/ 



where /i° x is the excess chemical potential of particles of 
diameter Cj. The first equality of Eq. ^ is related to the 



2 



reversible work needed to create a cavity large enough to 
accommodate a particle of infinite diameter Q. Note 
that the second equalities in Eqs. ([2]) and ((3]) are valid 
only if — > 3>({6*})- Internal consistency between 
Eqs. ([2]) and [3] implies 

(1 - Ss) _ = £o -* + g £ „_. ,4) 

Equation ((4| is independent of the scaling relation ([I]). 
Inserting the latter into the former one obtains the 
linear partial differential equation H(£s,z) = (1 — 
6)0£ 3 -ff (6i z)—zd z H(£s,z), whose solution is H(£ 3 , z) = 
-\E' I 1 * ) , where the function ^(y) remains undeter- 



mined. 

Therefore, if the (bulk) excess free-energy density of 
the HS mixture depends on the partial densities only 
through the four SPT variables, then it must necessarily 
have the simple scaling form 



$({6}) = -6 Mi -6) 



66 
i-6 



where we have called 



y 



CI 



127rCi(l-6) 
The pressure follows from either Eq. or (J3J) as 

6 66 



(5) 



(6) 



Pp. 



1-6 



(i-6) : 



[*(j/) + j/*'(z/)], (7) 



where ^'(y) = d*S(y)/dy. Agreement with the ex- 
act second and third virial coefficients simply requires 
&(y) = 1 + \y + 0(y 2 ). The excess chemical potential of 
component i is obtained from Eq. (0 as 



|i = - ln(l - &) + - A- [*(y) - £ 

opi 1-6 2 



[*(y) + 2tf*'(y)]7ro? + 0p^. 



i-e 



The unknown function ^(y) can be univocally derived 
from the one-component free-energy density. Setting 
(Tj = a, Eq. ([5]) reduces to 



,(»?) = -ln(l-r,) + 3 



1—7/ \ 1 — T] 



n 



(9) 



where (f> s (tl) = is the (bulk) excess free energy per 

particle of the one-component system. From Eq. (|9|) one 
gets 



*G0 



1 

3y 



y 



l + y 



- Mi + v) 



(10) 



The combination of Eqs. ([5]), and (fTU)) constitutes 
the main result of this Rapid Communication. It shows 



that, given the free energy of the pure HS system, the free 
energy of the HS mixture (if assumed to be truncatable 
with K = 3) is unambiguously known. In a more explicit 
form, 



*({U) 



6Mi-6)+4tt 
y 



6 2 



l + y 



6 

Mi + v) 



Analogously, Eq. (JT)) yields 

6 , „_ 6 2 



1-6 



47T 



6(1-6) 



l+y 



l + y 



(11) 



(12) 



where Z s {j]) is the compressibility factor pp/p of the 
one-component system. According to Eqs. (|lll) and 
(|12p. the thermodynamic properties of a HS mixture 
of total packing fraction n can be expressed in terms 
of those of a pure HS fluid with an effective packing 
fraction Vct{ = y/(l + y) = [1 + 12tt(1 - 6)6Al] _1 = 
[l + (l-r/)MiM 3 /r7M 2 2 ]- 1 < r). This contrasts with other 
approaches (2~H in which the reference one-component 
fluid has the same packing fraction (r) e g = if) as the mix- 
ture. 

It is worth noting that, while Eq. (|T2"j) repro- 
duces the exact second and third virial coefficients, 
it predicts a fourth virial coefficient given by B4 = 
(tt/6) 3 [Ml + 9M 3 Mf + 9MiM 2 Mf + (64 - 19)Mf/Mi] , 
where 64 is the reduced fourth virial coefficient of the 
one-component fluid. This expression is not consistent 
with a polynomial dependence on the mole fractions, 
B± = J2i,j.k.i x * x J x k x eBijke, except in the SPT case 
(64 = 19). This is a consequence of the approximate 
character of the FMT (or truncatability) ansatz at the 
level of the fourth virial coefficient, as already pointed 
out by Blaak [25|. 

A milder proposal consists of making the formal change 
*(y) -> A (6) + ^1(6)2/ in E q- ©, where the func- 
tions .Ao(6) an d ^4i(6) satisfy the differential equation 
(1 - 6)^0(6) + 6^1(6) = with the initial condi- 
tions A o (0) = 1, .Ai(O) = \. Although, in the general 
case, this breaks the consistency condition (QJ, the latter 
is verified in the one-component limit ^|/127r^i^3 — > 1. 
Given a desired compressibility factor Z s (r\)^ the differ- 
ential equation is closed with the algebraic relation (1 — 

6M0 (6) +26^1(6) = (i-6) 2 [(i-6)^(6)-i]/36- 

Choosing for Z a {rj) the CS EOS yields the WB mark II 
free energy, which was proposed in Refs. [HI 03 by a 
different method. 

The inhomogeneous case. Let us start by writing 
*(jv) = 1 + \y + A(y), where A(y) = 0(y 2 ), so that 
Eq. ([5]) becomes 



$ = $1 + $2 + *3 + $2A 



2£j 



(13) 



where $1 = -6 M 1 " 6), $2 =66/(1- 6), and $ 3 = 
^2 /247r(l ~ 6j) 2 ■ The SPT free energy corresponds to the 
linear approximation A(y) = 0. 
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In Rosenfeld's original FMT [l[ , the excess free-energy 
functional $[r, {^(r')}] in inhomogeneous situations can 
be constructed from the SPT bulk quantity [Eq. (fl"3|) 
with A(y) = 0] by following two basic steps. First, the 
four bulk SPT variables {6161676} are replaced by 
the weighted densities 



n a (r) = J dr ' Pi( 



r - r )oj n 



i(r'), 



(14) 



where the (scalar) weight functions are wo,j(r) = 
uj 2 ,i(r)/ircrf, wi,i(r) = u 2 s{v) /2-Kai, u) 2 s(v) = 8{r-a l /2), 
and W3,j(r) = Q(ai/2 — r). Second, to those four 
scalar weighted densities two vector weighted densities 
ni(r) and 112(1") are added. They are defined simi- 
larly to Eq. (HU, but with the vector weight functions 
wi j j(r) = oji.i(r)r/r and <jj 2 ^{r) = o> 2i i(r)r/r. The two 
vector densities ni and ri2 vanish in the bulk and so they 
are absent in the bulk free energy. A simple way of in- 
cluding ni and ri2 in 3>(r) and preserving the exact low- 
density behavior consists of making the formal changes 
66 ->■ nin 2 - ni • n 2 and £f -> n 2 (n\ - 3n 2 • n 2 ) in $ 2 
and $3, respectively. Therefore, Rosenfeld's functional is 
given by the right-hand side of Eq. (|T3"|) with A(y) = 
and 



*i({n a (r)}) = -n (r) ln[l - n 3 (r)], 



(15) 



* 2 (K(r)}) = "i(r)n a (r)-ni(r)-n 3 (,) > 



* 3 ({n«(r)}) 



1 - n 3 (r) 



n 2 (r)[n 2 (r) - 3n 2 (r) ■ n 2 (r)] 



24tt[1 - n 3 (r)] 



(17) 



In order to recover the exact one-dimensional (ID) 
functional for one-component systems ( "dimensional 
crossover"), Tarazona [4[ proposed to modify Eq. p7|) 

as 



<t>3(M) = 



n 2 {n\ 



3n 2 -n 2 ) + |[n 2 -n 2 -n 2 -Tr(n3)] 



24tt(1 - n 3 ) 



(18) 

where the additional tensor weighted density n 2 (r) is de- 
fined by Eq. (fT4]) with the tensor weight function 0, Q 
a?2,«(r) = w 2,i(r) (rr/r 2 — 1/3), I being the unit tensor. 
In the special case of a one-component system (cr, = er) 
confined to a ID geometry, i.e., Pi(r) = Xip^ 1D ^S(x)S(y), 
one has n 2 = p( 1D >er6(cr/2 - t)/ v / (cr/2) 2 - t 2 , n 2 /n 2 = 
(2/cr)t, and n 2 /n 2 = (4/cr 2 ) (tt — t 2 zz) +zz— |l, where 
t = \/.e 2 + y 2 and t = ix + yy. Insertion into Eq. f| 18[) 
yields $3=0 and spatial integration over the remaining 
two terms $1 + <I> 2 provides the exact ID free energy Q. 

In the case of the general class of consistent bulk free- 
energy densities (|13[) , it seems quite natural to construct 
the corresponding class of functionals by applying the 
replacements (|15|) and (|T5)) , together with either Eq. (fTT)) 
(vector densities) or Eq. (fTg]) (tensor densities) . In order 



to specify a particular functional one only needs to choose 
a thermodynamic description for the pure fluid as input 
and then obtain the function ty(y) = l + \y+A{y) via Eq. 
© or, equivalently, Eq. (JTUJ). Since A(y) = 0(y 2 ) and 
$3, as given by Eq. (fl"8|) . vanishes in the one-component 
ID limit, it is obvious that the whole class of free-energy 
densities (fT3|) becomes exact in that limit. On the other 
hand, the presence of the nonlinear term A(y) in Eq. 
(|13[) induces in the associated direct correlation function 
a spurious simple pole at r = [26j . 

Application to bulk properties. As an application of 
our scheme, let us assume that the bulk excess free energy 
per particle of the pure fluid has the form 



4>s(v) = -ln(l -rj) + 3n 
1 



1 - r,/2 



+A 



v- 



(i-v) 2 



(1-'/) 

The associated EOS of the pure HS fluid is 



z s (v) = 



1 + r] + r) 2 - A?/ 3 



(19) 



(20) 



The forms (fT^|) and (j2"0")) encompass the virial-route 
Percus-Yevick (PY) (A = 3), the compressibility- route 
PY (or SPT) (A = 0), and the CS [H (A = 1) EOSs. In 
general, the coefficient A is related to the reduced fourth 
virial coefficient by 64 = 19 — A. Insertion of Eq. p^|) 
into Eq. (JTOj) yields 



1 A 
1 + -y+ - 
2 y 3 



1 In 1 + y) 

2 y 



(21) 



In the SPT case (A = 0), ^{y) is simply approximated 
by the exact expansion in powers of y truncated after 
the linear term. In general, Eq. (f2"Tj) gives ^(y) = 1 + 
hy — (A/9)y 2 + 0(y 3 ). According to Eq. ([5]), the unique 
frcc-cncrgy density of the mixture consistent with the 
one-component expression (|19p is 



A\ 66 



<&({6}) = -6Mi-6) + 1 + 7 



1-6 



1 - 



A 



3 j 24tt(1 -i z ) 2 

«1„ 1 , ?2 



-4ttA^- In 

6 



1 



12tt6(1-6) 



(22) 



The associated EOS can be obtained from either Eq. ([2]) 
or Eq. ©. The result is 



6 



1-6 



A 

1+ 3 



66 



(1-6) 2 



3; 127r(l-60 3 

4ttA 



6 2 6 



6 ^ + 12tt6(i-6 



(23) 
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It is easy to check that Eqs. ([22]) and ([23)1 reduce to Eqs. 
(Unj) and ([2"0"]) , respectively, in the one-component limit 
cr t cr. 

If we take as input the SPT thermodynamic descrip- 
tion for the one-component fluid, i.e., Eq. (jT9|) with 
A = 0, then the output free energy for the mixture given 
by Eq. ([2"2"]) coincides with the true SPT result for mix- 
tures. This is an expected result since, as is well known, 
the SPT is thermodynamically consistent with Eqs. ([2]) 
and ©. On the other hand, the SPT EOS is known to 
overestimate the pressure of the HS fluid. A much bet- 
ter description is provided by the CS description, i.e., by 
Eqs. (HU) and ([20]) with A = 1. Taking this as input, 
the consistent extension to mixtures is provided by Eqs. 

and (22]) with A = 1. The result differs from the 
most popular extension of the CS equation to mixtures, 
namely the BMCSL theory 0Q. The BMCSL EOS 
is 



Pp. 



6 



66 



1-6 (1-6) 



6 3 (3-6) 

1 367r(l-£ 3 ) 3 ' 



(24) 



According to the thermodynamic relation ([2]), the free- 
energy density corresponding to Eq. ([24]) is 



$({6}) = -6in(i-6) 



6 3 



36< 3 2 



6 



(1-6) 



66 
1-6 

' Mi - 6) 



(25) 



On the other hand, if Eq. ([3]) is used instead, the result 
is 



-6M1-6 ) + -, r + ^F71 Txi' 26 

1-6 36tt(1-6) 



The difference between Eqs. ([25]) and ([26]) reflects that, 
in contrast to Eq. ([23]). the BMCSL EOS ([24]) is inconsis- 
tent with respect to the simultaneous verification of the 
thermodynamic relations © and ([3]). The free-energy 
density ([25]) was the one used in Ref. [f| to construct the 
WB version of the FMT. 

The alternate EOS {23), as well as the BMCSL, PY, 
and SPT EOSs, share the property that, at a given to- 
tal packing fraction, the compressibility factor /3p/p de- 
pends on the size distribution only through the dimen- 
sionless combinations of moments 71 = M1M2/M3 and 
72 = M2/M3. On the other hand, while the dependence 
of the PY, SPT, and BMCSL EOSs on both 71 and 72 
is linear, the last term on the right-hand side of Eq. ([2"3"]) 
introduces a nonlinear dependence. 

As a test of the superiority of this free-energy density, 
Eq. ^ with A = 1, over the BMCSL one, Eq. ([23]). 
Table [I] compares the respective deviations of the theo- 
retical values with respect to the simulation ones [27j for 
the compressibility factor f3p/p and the (reduced) chem- 
ical potentials /3pi and /3/i2- The system chosen is the 
one with the highest packing fraction (7/ = 0.49) and the 
largest size disparity (o^/ci = 0.3) considered in Ref. 



TABLE I. Difference between theoretical and simulation val- 
ues [27| for the compressibility factor pp/p and the (reduced) 
chemical potentials /3pi and fipi in a binary mixture with a 
total packing fraction r\ — 0.49, a size ratio o"2/a"i = 0.3, and 
several mole fractions xi. 



HPp/p) 



Eq. (22) Eq. (25) 


Eq. (22) Eq. (25) 


Eq. (22) 


Eq. (25) 


-0.02 


-0.08 


-0.1 -0.8 


-0.01 


-0.04 


-0.01 


-0.07 


-0.05 -0.4 


-0.002 


-0.03 


-0.007 


-0.05 


-0.0002 -0.2 


0.001 


-0.02 


-0.01 


-0.04 


-0.006 -0.07 


0.004 


-0.01 


-0.02 


-0.04 


-0.03 -0.05 


0.005 


-0.01 



[27j . We observe that the deviations are strongly reduced 
by the use of the consistent free energy (f2"2"[) . Compari- 
son with Table I of Ref. [l7| shows that the WB mark II 
free energy predicts even better values for the pressure. 
However, Eq. ([22]) predicts more accurate values of fXi 

For the largest 



111 

8 ' 4 ' 2 ' 



for x\ = j and of p,i for x\ 

mole fraction {x\ = |), Eq. ([22]) and the WB mark II free 
energy give practically the same results for p, pi , and /i 2 ■ 

It is interesting to remark that Eq. ([22]) with A = 1 
does predict a (metastable) demixing transition, in con- 
trast to the SPT [Eq. ^ with A = 0] and the BM- 
CSL [Eq. ([25]) ] descriptions. In particular, the predicted 
coordinates of the critical point are (x\ c , n c , /3p c a{) = 
(0.0022,0.686,7305) and (0.000 32,0.594,19 404) for the 
size ratios (72 Jo 1 = and ^ , respectively. 

Conclusions. To sum up, once the FMT (or trun- 
catability) ansatz for the bulk free-energy density is as- 
sumed, basic consistency conditions impose that it must 
necessarily have the functional form expressed by Eq. 
([5]); where the scaled variable y is given by Eq. (6]) and 
the dimensionless scaling function ^(y) remains arbi- 
trary, except for the low-density requirement ^{y) = 
1 + \y + 0(y 2 ). This arbit rarincss, however, can be 
exploited to seek consistency with any desired EOS in 
the one-component case via Eqs. © or ([TU]) . Thus, Eq. 
([5]) represents a class of consistent FMT free energies, 
parametrized by the one-component function. This class 
includes the SPT free energy [Eq. ([22]) with A = 0] as a 
particular case. However, the WB free energy ([25]) docs 
not belong to the class. In fact, the unique member of the 
class consistent with the CS EOS [Eq. (g2]) with A = 1] 
turns out to be much more accurate than the WB free 
energy, as shown by Table U Of course, other choices 
are possible and it can be reasonably expected that the 
more accurate the one-component input, the better the 
multicomponent output. 

While most of the Rapid Communication has dealt 
with the bulk free energy, the extension of the re- 
sults to confined geometries and/or to external poten- 
tials is straightforwardly given by Eq. (]13[) complemented 
by Eqs. ([IS]) and ([16]) plus Eq. (TTT]) (including vector 
weighted densities) or Eq. ([18]) (including tensor weighted 
densities). 
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As any other approximation, the proposal presented 
here has strengths and weaknesses. Weak points are 
that, except for the genuine SPT free energy, the virial 
coefficients beyond the third one do not depend polyno- 
mially on the mole fractions and the direct correlation 
function stemming from the inhomogeneous functional 
possesses a pole at r = 0. These shortcomings arc reflec- 
tions of the approximate character of the FMT ansatz 
®({Pi}) On the other hand, strong features 

are the internal consistency with Eqs. ([1]) and (|1]), the 
flexibility to accommodate any desired one-component 



thermodynamic description, and the accuracy of the bulk 
free energy when an accurate one-component theory is 
chosen as input. Moreover, preliminary calculations for 
confined fluids provide encouraging results [28| . 
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